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Abstract 

The potential benefits of multiple-antenna systems may be limited by two types of channel degrada- 
tions — rank deficiency and spatial fading correlation of the channel. In this paper, we assess the effects 
of these degradations on the diversity performance of multiple-input multiple-output (MIMO) systems, 
with an emphasis on orthogonal space-time block codes, in terms of the symbol error probability, 
the effective fading figure (EFF), and the capacity at low signal-to-noise ratio (SNR). In particular, we 
consider a general family of MIMO channels known as double -scattering channels, which encompasses a 
variety of propagation environments from independent and identically distributed Rayleigh to degenerate 
keyhole or pinhole cases by embracing both rank-deficient and spatial correlation effects. It is shown that 
a MIMO system with rix transmit and riR receive antennas achieves the diversity of order max(rtT^ws"^TtR) 
in a double-scattering channel with ns effective scatterers. We also quantify the combined effect of the 
spatial correlation and the lack of scattering richness on the EFF and the low-SNR capacity in terms 
of the correlation figures of transmit, receive, and scatterer correlation matrices. We further show the 
monotonicity properties of these performance measures with respect to the strength of spatial correlation, 
characterized by the eigenvalue majorization relations of the correlation matrices. 

Index Terms 

Channel capacity, diversity, double scattering, fading figure, keyhole, multiple-input multiple-output 
(MIMO) system, orthogonal space-time block code (OSTBC), spatial fading correlation, symbol error 
probabiUty (SEP). 

I. Introduction 

Recent rapid advances in multiple-input multiple-output (MIMO) communication theory and 
growing cognizance of the tremendous performance gains achieved by MIMO techniques [l]-[9] 
have spurred efforts to integrate this technology into future wireless systems such as wireless local 
area networks (WLANs) and 4G cellular systems. One of the approaches to exploiting diversity 
capability of MIMO chaimels is the use of orthogonal space-time block codes (OSTBCs), 
which have drawn considerable attention because they attain full diversity with scalar maximum- 
likelihood (ML) decoding [7]-[9].i 

'However, OSTBCs with arbitrary complex constellation cannot provide tiie full diversity and full transmission rate 
simultaneously for more than two transmit antennas [8, Theorem 5.4.2] (see also [10]-[13]). A new class of quasi-orthogonal 
codes has been proposed in [14]-[16] with the tradeoff between the decoding complexity, transmission rate and/or diversity. 
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In general, the potential benefits of multiple- antenna systems may be limited by rank deficiency 
of the channel due to double scattering or the keyhole effect, for example, as well as spatial fading 
correlation due, for instance, to insufficient spacing between antenna elements [17]-[30]. Some 
mechanism rendering a MIMO channel rank deficient cannot be explained by the archetypal 
model based on single-scattering processes [26], [27]. To address this issue, a double-scattering 
MIMO model has been proposed recently in [24] wherein the channel matrix is characterized by a 
product of two statistically independent complex Gaussian matrices, in contrast to the common 
single complex Gaussian matrix characterization for wireless MIMO channels.^ This double- 
scattering model can capture both rank-deficient and spatial correlation effects of MIMO channels 
and encompass a variety of propagation environments, bridging the gap between an independent 
and identically distributed (i.i.d.) Rayleigh case and a degenerate one-rank channel known as a 
keyhole or pinhole channel. There are other recent attempts to modeling MIMO channels for 
more realistic scattering environments (e.g., double or multibounce diffuse scattering) beyond 
single scattering [31]-[34]. 

The effects of rank deficiency and spatial correlation on the capacity of MIMO channels 
are relatively well understood (see, e.g., [17]-[30]). From a capacity point of view, it has been 
known that at high signal-to-noise ratio (SNR), the spatial fading correlation reduces the diversity 
advantage — a parallel shift of the capacity curve over SNR in decibels (dB) — offered by multiple 
antennas, whereas the rank deficiency decreases the spatial multiplexing benefit — a slope of the 
capacity curve over SNR — of multiple- antenna channels [21]. Previously, the performance of 
space-time coding in the presence of spatial fading correlation has been extensively studied for 
the most popular Rayleigh, Rician, and Nakagami-m fading [35]-[40]. Also, the effect of rank 
deficiency has been investigated in [41]-[44] for a special case of the keyhole channel. 

The objective of this paper is to assess the effects of double scattering on the diversity 
performance of MIMO systmes in a communication link with n-r transmit antennas, tir receive 
antennas, and ns effective scatterers on each of the transmit and receive sides, which is referred 
to as a "double-scattering (rix, ns, nR)-MIMO channel." Due to the channel decoupling property, 
the OSTBC converts a MIMO fading channel into identical single-input single-output (SISO) 
subchannels, each for a different transmitted symbol, with a path gain given by the Frobenius 

^In [24], the model was validated by simulations using ray tracing techniques. 
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norm^ of the channel matrix H [38]-[42]. As a resuh, the maximum achievable diversity perfor- 
mance of MIMO systems can be characterized by the statistical property of ||-ff||p. Therefore, 
using the OSTBC as a pivotal MIMO diversity technique"^ (particularly, in the absence of channel 
knowledge at the transmitter), we analyze the relevant performance measures in double- scattering 
(nx, ns, nR)-MIMO channels, namely: i) the symbol error probability (SEP) [49], ii) the effective 
fading figure (EPF) [50]-[52], and iii) the capacity in a low-SNR regime [53], [54]. 

Diversity in communication can ameliorate system performances in behalf of error probability, 
information rate, and signal fluctuation due to fading. Prom a error probability viewpoint, the 
diversity attacks a high-SNR slope of the SEP curve, i.e., diversity order. In contrast, the diversity 
(from a capacity point of view) affects a low-SNR slope of the capacity curve rather than a high- 
SNR slope. Por example, the high- and low-SNR slopes (bits/s/Hz per 3 dB) of the capacity for 
i.i.d. Rayleigh-fading MIMO channels are given by 



respectively [53]. While the high-SNR capacity slope Soo is limited by the spatial multiplexing 
gain min (nx, ?^r), the low-SNR capacity slope 5*0 is limited by the diversity gain amounting to 
the harmonic mean of and n^. Therefore, the capacity is multiplexing-limited in the high-SNR 
regime, but is diversity-limited in the low-SNR regime. At high SNR, the diversity advantage 
serves only to provide the power offset (i.e., the parallel shift of the capacity curve) [21]. These 
lessons stimulate a shift of focus to the low-SNR regime in analyzing the diversity effect on the 
capacity behavior. More inherently, diversity systems aim to reduce signal fluctuations due to 
the nature of fading. The EPP measure is defined as a variance-to-mean-square ratio (VMSR) 
of the instantaneous SNR (see Definition [T]). This quantity can be used to assess the severity 

'The Frobenius norm of an m x n matrix A = (Atj) is defined as 



where tr(-) and f denote the trace operator and the transpose conjugate of a matrix, respectively. 

''if the transmitter has channel knowledge, the maximum MIMO diversity can be achieved by transmit beamforming (often 
called maximum ratio transmission (MRT) or MIMO maximal-ratio combining) in the eigenspace of the largest eigenvalue of 
the Gramian matrix H'' H [45]-[48]. 



5*00 = min(nT,nR) 



2nTnR 



— \ 
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of fading and the effectiveness of diversity systems on reducing signal fluctuations. The main 
results of this paper can be summarized as follows. 

• We show that the achievable diversity is of order 

max(nT,ns,nR) ' 

Hence, if the channel is "rich-enough," that is, the number of effective scatterers is greater 
than or equal to the numbers of transmit and receive antennas, the full spatial diversity 
order of riTn^, can be achieved even in the presence of double scattering. 

• We derive exact analytical expressions for the SEP in three cases of particular interest: 

1) spatially uncorrected double scattering (includes i.i.d. and keyhole channels as special 
cases); 

2) doubly correlated double scattering (includes a spatially correlated MIMO channel 
where spatial correlation is present at both the transmitter and the receiver); 

3) multiple-input single-output (MISO) double scattering (corresponds to a pure transmit 
diversity system wherein a burden of diversity reception at the receive terminal is 
moved to the transmitter — original motivation of space-time coding [6]-[8]). 

• We derive the EPF and the low-SNR capacity of double-scattering (nx, ns, nR)-MIMO 
channels. The results show that these performance measures are completely characterized 
by the correlation figures of transmit, receive, and scatterer correlation matrices.^ 

• The EPP as a functional of the eigenvalues of correlation matrices is monotonically in- 
creasing in a sense of Schur (MIS).^ We show that the maximum possible increase in the 
EPP due to double scattering is a sum of correlation figures of the transmit and receive 
correlation matrices, which eventuates when the scatterers tend to be fully correlated or the 
keyhole propagation takes place, that is, when only a single degree of freedom is available 
in the channel for communications. 

• The low-SNR capacity slope as a functional of the eigenvalues of correlation matrices is 
monotonically decreasing in a sense of Schur (MDS). We also obtain the low-SNR capacity 
of a double-scattering MIMO channel without the constraint of orthogonal input signaling. 

'The correlation figure is defined as a ratio of the second-order statistic of the spectra of correlation matrices to that of the 
fully correlated matrix (see Definition |3- 

*See Appendix U for the notions of Schur monotonicity and majorization. 



DRAFT 



July 5, 2006 



SHIN AND WIN: MIMO DIVERSITY IN THE PRESENCE OF DOUBLE SCATTERING 



5 



This enables us to assess the penalty of the use of OSTBCs (for achieving full diversity 
with simple decoding) on spectral efficiency in the low-SNR regime. 

We note in passing that all the mathematical and statistical results (on the monotonicity in a 
sense of Schur and random matrices) obtained in the appendices are applicable to many other 
problems related to multiple-antenna communications — for example, capacity analysis of MIMO 
relay channels [5] and spatially correlated MIMO channels [21]-[23], and error probability 
analysis of multiple- antenna systems with cochannel interference [55], [56]. 

This paper is organized as follows. In Section im the system model considered in the paper is 
presented. Section [nil analyzes the achievable diversity and the SEP in the presence of double 
scattering. Section |lVl analyzes the EPF and the low-SNR capacity (with and without the use 
of OSTBCs) of double- scattering (nx, ng, r;,R)-MIMO channels. Section fVl concludes the paper. 
Apropos of our study, the notions of majorization and Schur monotonicity are briefly discussed in 
Appendix U In Appendix HH we provide supplementary useful results on some statistics derived 
from complex Gaussian matrices. 

Notation: Throughout the paper, we shall use the following notation. N, M, and C denote the 
natural numbers and the fields of real and complex numbers, respectively. The superscripts *, 
T, and f stand for the complex conjugate, transpose, and transpose conjugate, respectively. /„ 
and Omxn represent the n x n identity matrix and the m x n all-zero matrix, respectively. (Aij) 
denotes the matrix with the (i, j)th entry Aij and deti<jj<n (Aij) is the determinant of the nxn 
matrix (Aij). tr(A), etr (^) = e*^*^"^^, and \\A\\p denote the trace, exponential of the trace, and 
Frobenius norm of the matrix A, respectively. ® and © denote the Kronecker (direct) product 
and direct sum of matrices and vec (A) denotes the vector formed by stacking all the columns of 
A into a column vector. Also, we denote Ai^A2®---^Anhy <S>7=i -^i ^^'^ © A2 © ■ ■ ■ © A„ 
by With a slight abuse of notation, a positive-semidefinite matrix A is denoted by 

A > and a positive-definite matrix A is denoted hy A > 0. Finally, for a Hermitian matrix 
A E C"^" with the eigenvalues Ai, A2, . . • , A„ in any order, g (A) denotes the number of distinct 
eigenvalues of A. Also, A^^^ and (A), k = 1,2, . . . , g (A), denote the distinct eigenvalues of 
A in decreasing order and its multiplicity, respectively, that is, A^^^^ > A^2) > • • • > '^(^(A)) ^^'^ 
Y.iil r, {A) = n. 
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II. System Model 

We consider a MIMO wireless communication system with ut transmit and receive 
antennas, where the channel remains constant for an integer multiple of Nc (> n^) symbol 
periods and changes independently to a new value for each coherence time. We assume that the 
channel is perfectly known at the receiver but unknown at the transmitter. 

A. Orthogonal Space-Time Block Codes 

A space-time block coded MEMO system in double-scattering channels is illustrated in Fig. [T] 
During an Ai'c-symbol interval, symbols Xi E S, i = 1,2, . . . , N, are encoded by an OSTBC 
defined by an Nc x transmission matrix G, where S is two-dimensional signaling constellation 
[8], [9]. A general construction of complex OSTBCs with the minimal delay and maximal 
achievable rate was presented in [10, Proposition 2]. This construction of the OSTBC for n-p 
transmit antennas gives the maximal achievable rate [10, Theorem 1] 



n 



\\0g2 Ut] + 1 
2 [logs "tI 



(1) 



where [x] denotes the smallest integer greater than or equal to x. For example, Alamouti's code 
[ -x^ xl] is a one-rate OSTBC employing two transmit antennas [7] and 



Xi 



X2 





Xs 



n 





-Xs 
Xl 



(2) 



is a 3/4-rate OSTBC for four transmit antennas [10]. 



B. Signal and Channel Models 

For a frequency-flat block-fading channel, the tir x Nc received signal can be expressed in 
matrix notation as 

Y = Hg^ + W (3) 

where H E C""^""^ is the random channel matrix whose (i,j)th entries Hij, i = 1,2, . . . ,nR, 
j = 1,2, . . . , ut, are complex propagation coefficients between the jth transmit antenna and the 
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ith receive antenna with E{|/Jjjp} = 1, and W ~ ■AfnB,N^ i'^nB^xN^^ NqI n^, I is the complex 
additive white Gaussian noise (AWGN) matrix (see [21, Definition II. 1] and [21, (1)] for the 
definition and distribution of complex Gaussian matrices)/ The total power transmitted through 
nx antennas is assumed to be V and hence, the average SNR per receive antenna is equal to 

7 = V/No. 

For double- scattering {n^, no,, nR)-MIMO channels (see Fig. [T]), the channel matrix H can be 
written as [21], [24] 

F = ^$fF,$f/r,$f (4) 

where no, is the number of effective scatterers on each of the transmit and receive sides. Hi and 
H2 are statistically independent. Hi ~ A/'„j,,ns (0„j^xns, -^riR, -^ns)' -^2 ~ -Mig.nT (Ong 
and Hermitian positive-definite matrices ^s? and $r are rix x rix transmit, ns x scatterer, 
and riR x receive correlation matrices with all diagonal entries 1, respectively.^ This model 
can include the rank-deficient effect of MIMO channels as well as spatial fading correlation 
by controlling and the correlation matrices and Therefore, Q is a general 

family of MIMO channels spanning from the i.i.d. Rayleigh case (ns oo with $t = Iut^ 
= Ins^ = Inj^) to the degenerate keyhole or pinhole case (ns = 1 with $t = In-r, 
= Itib) [24]. Note that the separability of correlation in (|4l) is a generalization of the well- 
known 'Kronecker model' [17], [18]. Although there are some attempts to reporting discrepancy 
between this separable correlation model and physical measurements (see, e.g., [57], [58]), the 
Kronecker correlation model has been accepted widely due to its experimental validation from 
European Project [19] and analytical tractability. 

In [20], so-called stochastic rank deficiency — meaning that the channel is rank deficient 
due to fading correlation, i.e., the correlation matrices have zero eigenvalues — was deemed 
as an important feature when dealing with fading correlation. However, this form of channel 
degeneracy cannot cover the case where the channel exhibits rank deficiency even when fading 
is uncorrelated. In contrast, we shall restrict and $r to positive-definite (i.e., full rank) 

matrices in the paper. This implies that the rank of H is equal to min (nx, ns, n^) with probability 
one. Therefore, rank deficiency can be distinguished from the fading correlation effect and may 

'There exist minor typos in [21, Definition II. 1]; the covariance matrix S ® * should be read as S"^ (g) 9. 
^In general, a correlation matrix is positive semidefinite with all diagonal entries 1. 
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occur only due to the lack of scattering richness with ns less than min(nT,nR). This also 
enables us to discriminate a one-rank /aZZ); correlated scenario from a degenerate keyhole MIMO 
channel [29], and grants the channel to exhibit rank deficiency with uncorrelated fading (e.g., 
ns < min (nx, n^,) with $t = 1^, = ^ng, $r = Inn)- 
Let Hi = and E2 = ^g^^Ha^T^^ then we have 

H=^E,E, (5) 



where Hi ~ J^nn,ns (0„Rxns,*R,^ns) and E2 ~ A^„s,„t (0„sxnT'*s,^T) are statistically inde- 
pendent complex Gaussian matrices. 

III. Symbol Error Probability 

With perfect channel knowledge at the receiver, orthogonal space-time block encoding and 
decoding convert a MIMO fading channel into equivalent SISO subchannels, each for a 
different symbol, with a path gain \\H\\p [38]-[42] (as shown in Fig. [T]). Consequently, the 
performance of OSTBCs is completely characterized by the statistical behavior of ||^||p and the 
instantaneous SNR for each of the SISO subchannels, denoted by 7stbc5 is given by [41], [42] 

l\\H\\l 

TSTBC = (O) 

nT/v 

To evaluate the SEP, we need the probability density function (pdf) or the moment generating 
function (MGF) of 7stbc- For double- scattering (nx, ns, nR)-MIMO channels, the MGF of 7stbc 
can be written as 

0,,,3, (s;7)=E|etr (^--^HH^ 
Eh, h, < etr , ^ 




(7) 
(8) 

where Q and (EJ) follow from Lemma [T] in Appendix HH 
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A. Achievable Diversity 

Before devoting to deriving the SEP expressions, we discuss the diversity order achieved by 
the OSTBC. In general, the achievable diversity order can be defined as 

d A lim (9) 

7^00 log 7 

where Pg denotes the SEP for two-dimensional signaling constellation with polygonal decision 
boundaries. In the absence of double scattering, the OSTBC provides the maximum achievable 
diversity order of rix^R- The corresponding diversity order in double- scattering (nT,ns,nR)- 
MIMO channels is given by the following result. 

Theorem 1: The diversity order achieved by the OSTBC over double- scattering {n'T,ns,n^)- 
MIMO channels is 

OsTBC = 7 r- (10) 

max(nT,ns,nR) 

Proof: See Appendix IIII-Al □ 

Theorem [l] states that if the number of effective scatterers is greater than or equal to the 
numbers of transmit and receive antennas, the OSTBC provides the full diversity order of utUb^ 
even in the presence of double scattering. 

We now present analytical expressions for the SEP of the OSTBC for three cases of particular 
interest — spatially uncorrelated double scattering, doubly correlated double scattering, and MISO 
double scattering. In what follows, a spatial correlation environment of double-scattering channels 
is denoted by T = ($t,^S5^r) for given nx, ns, and n^. 

B. Spatially Uncorrelated Double Scattering 

Consider a spatial correlation environment Tuc = (/„^ , /„g , /„j, ) . This spatially uncorrelated 
double-scattering scenario includes i.i.d. and keyhole MIMO channels as special cases. 
Let ni = min (rix, ^s), ^2 = max (nx, ris), and the rii x rii random matrix T be 

E2EI, if ns < ut 

t (11) 
^2^2, if ns>nT, 

which is a matrix quadratic form in complex Gaussian matrices [21, Definition n.3]. Then, from 
© and (fT47t in AppendixInD the SEP of the OSTBC with M-PSK signaling in double-scattering 
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(nx, ^s, ?t,r)-MIMO channels can be readily written as 



1 r^l.Jr , 91 



e.MPSK 



, E<^det + ^^-^$j,®T }d9 (12) 




where we have used the fact that H2H2 and H2H2 have the same nonzero eigenvalues.^ 

In the absence of spatial correlation, the random matrix T has the Wishart distribution 
VV„^ (?22,/ni) [21, Definition 11.2]. Applying Corollary El in Appendix HH to ([Hb . we obtain 
the SEP for this spatially uncorrelated environment Tuc as 

^"^MPSK = ^ [ det {G-''^ {6)] de (13) 

where 

ni 

A"^-'^' = W{n2-k)\{k-l)\ (14) 

k=l 

and G"'^"'*'* {6) = (G"j"'*** (6)) is the rii x rii Hankel matrix whose («, j)th entry is given by 

G^;-'^^ (9) = (n2 - ni + ^ + J - 2)! 2^0 (n2-ni + t + j- 1, n^; ^^—TT^l ■ d^) 

\ usUt /V sm u J 

Example 1 (Uncorrelated Extremes — Keyhole and I.I. D.): The i.i.d. and keyhole MIMO chan- 
nels are two extreme cases of spatially uncorrelated double scattering (i.e., ns = 00 and ns = I, 
respectively). If ns = 1, then ni = 1 and ^2 = ut- Hence, (fT3t reduces to [41, eq. (11)] for 
keyhole MIMO channels. As 00, (fT3t becomes [42, eq. (26)] (with a Nakagami parameter 

m = 1) for i.i.d. Rayleigh-fading MIMO channels. 

Fig. El shows the SEP of 8-PSK (2.25 bits/s/Hz) versus the SNR 7 in spatially uncorrelated 
double-scattering (4, ns, 2)-MIM0 channels when ns varies from 1 (keyhole) to infinity (i.i.d. 
Rayleigh). We can see that as ns increases, the SEP approaches that of i.i.d. Rayleigh-fading 
MIMO channels in the absence of double scattering. This resembles the behavior in Rayleigh- 
fading channels with diversity reception, that is, the channel behaves like an AWGN channel 

'As mentioned in the proof of Theorem Q The SEP for the general case of arbitrary two-dimensional signaling constellation 
with polygonal decision boundaries can be written as a convex combination of terms akin to <147> . Thus, our results can be 
easily extended to any two-dimensional signaling constellation. 
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(diversity order of oo) as the number of receive antennas increases. Observe that when no, > 4, 
the slope of the SEP curve at high SNR is identical to that of the i.i.d. case. This example 
confirms the result of Theorem [T] the diversity orders are equal to c^stbc = 2, 4, and 6 for 
ns = 1, 2, and 3, respectively, whereas (isTBc = 8 for = 5, 10, 20, 100, and oo (i.i.d.). 
A clearer understanding about the diversity behavior is obtained by referring to Fig. |3l where 
the SEPs of 16-PSK Alamouti (4 bits/s/Hz) and (3 bits/s/Hz) OSTBCs versus the SNR 7 
in spatially uncorrelated double- scattering (^t, ns, r;,R)-MIMO channels are shown. Using (fTOI) . 
we can easily show that the Alamouti and ^4 codes achieve the diversity order of c?stbc = 2 for 
(2, 3, 1) and (4, 2, 1) channels; c/stbc = 6 for (2, 5, 3) and (4, 3, 2) channels; and c^stbc = 20 for 
(2, 10, 11) and (4, 5, 5) channels. As can be seen, we obtain a close agreement in the slopes of 
the SEP curves, corresponding to the same value of c/stbc^ at high SNR. 

C. Doubly Correlated Double Scattering 

Consider a spatial correlation environment T^c = (^t, -^ng) ^r)' where spatial correlation 
exists only on the transmit and receive ends. Note that this scenario includes a spatially correlated 
MIMO channel in the absence of double scattering (ns = 00) as a special case. Let and A^, 
i = 1,2,..., TT-T, j = 1, 2, . . . , n^, be the eigenvalues of $x and $r in any order, respectively. 
Suppose that ris > ^t- Then, T ~ VVn^ (ns,$T)- Applying Theorem flOl in Appendix |ll| to (fT2t . 
we obtain the SEP in the environment T^c as 

^s'!M'psK = ^^^°det([Gr^(^) Qt''{0) ■■■ Q':(ilM)dO (16) 

with 

Where B^^ = {^^) and G^^ (0) = (G^,-f (0)), k = 1, 2, . . . , ($t), are x r, ($t) 
matrices whose (i, j)th entries are given respectively by 

Bgf = (-irM^-J + l),-iA5;-^+^ (18) 
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and 



(?(*r) •rp(*R) C 




G£f (^) = E E '^^.'^ (^i^) ■ ^ 



p=l q=l (. 



UT + i + j -2)1 



X 2F0 \ ns - riT + i + j - 1, q 



usUtTI sin 





(19) 



In ([T9l) . q ($r) is the {p, g)th characteristic coefficient of (see Definition |4| in Appendix Hll). 

Fig. m shows the SEP of 8-PSK versus the SNR 7 in doubly correlated double-scattering 
(4, 10, 4)-MIM0 channels. In this figure, the transmit and receive correlations follow the constant 
correlation $t = = ^4^^ {p), defined by (l53t in Appendix U and the correlation coefficient 
p ranges from (spatially uncorrelated double scattering) to 0.9. The characteristic coefficients 
of the constant correlation matrix are given by (11311) and (11321) (see Example |6l in Appendix HH). 
For comparison, we also plot the SEP of i.i.d. Rayleigh-fading MEMO channels. In Figure El we 
can see that the SNR penalty due to double scattering with ns = 10 (in the absence of spatial 
correlation) is about 1 dB at the SEP of 10^ and it becomes larger than 2.5 dB for p > 0.5. In 
Fig. 121 the SEP of 8-PSK ^4 at 7 = 15 dB is depicted as a function of a correlation coefficient 
p for doubly correlated double- scattering (4, ris, 4)-MIM0 channels with constant correlation 
$T = ^R = ^4''^ (p) when ns = 5, 10, 20, 50, 100, and 00 (doubly correlated Rayleigh). This 
figure demonstrates that double scattering and spatial correlation degrade the SEP performance 
considerably. 

D. MISO Double Scattering 

Finally, we consider a double- scattering MISO channel. This is a pure transmit diversity system 
wherein the burden of diversity reception at the receive terminal is moved to the transmitter. 
The SEP in double-scattering MISO channels can be obtained from (jH) with = 1 as 



Let Xf, i = 1, 2, . . . , ns, be the eigenvalues of $s in any order. Then, applying Theorem [TTI in 




(20) 
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Appendix ini to (l20b . we obtain 



1 





) 



PI 



jmiso-ds 

e.MPSK ~ 



TT 



E E E E '^p. ($s) A",, ($t) / ^, J 



nsn-rTZ sin^ 9 



dO (21) 



p=l q=l i=l j=l 



where Xp i ($s) and Xq j ($t) are the characteristic coefficients of $s and respectively. 

The effects of the spatial correlation and the number of effective scatterers on the SEP 
performance in MISO channels can be ascertained by referring to Fig. |^ where the SEP of 8- 
PSK ^4 at 7 = 25 dB versus ns is depicted for double- scattering (4, ns, 1)-MIM0 channels. The 
transmit and scatterer correlations follow the constant correlation $t = ^4^^ (p) and $s = ^^^'^^ (p) 
where p varies from to 0.9. Note that the maximum achievable diversity order is equal to 
o^STBC = 4 for ns > 4. Hence, the SEP performance improves rapidly as ns increases, and 
approaches the corresponding SEP in the absence of double scattering. 

IV. Effective Fading Figure and Low-SNR Capacity 

In this section, we access the combined effect of rank deficiency and spatial correlation on the 
performance of OSTBCs in terms of the EFF and the capacity in a low-SNR regime. It will be 
apparent that these performance measures are completely characterized by the kurtosis of ||p. 

A. Effective Fading Figure 

One of the goals of diversity systems is to reduce the signal fluctuation due to the stochastic 
nature of multipath fading. Therefore, it is of interest to characterize the variation of the instan- 
taneous SNR at the output where the amount of signal fluctuations is measured. The following 
measure can be used to assess the severity of fading and the effectiveness of diversity systems 
on reducing signal fluctuations. 

Definition 1 (Effective Fading Figure): For the instantaneous SNR 7 at the output of interest 
in a communication system subject to fading, the effective fading figure (EFF) in dB for the 
output SNR 7 is defined as the VMSR of 7, i.e.. 



It should be noted that the EFF is akin to the notions of the normalized standard deviation 
(NSD) of the instantaneous combiner output SNR [50]-[52] and the amount of fading (AF) 




(22) 
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[59], [60]. The AF, as defined in [59, eq. (2)], is purely to characterize the amount of random 
fluctuations in the channel itself and conveys no information about diversity systems. In contrast, 
the NSD is a measure of the signal fluctuations at the diversity combiner output, enabling us to 
compare the effectiveness of diversity combining techniques such as maximal-ratio combining 
(MRC), equal-gain combining (EGC), selection combining (SC), and hybrid section/maximal- 
ratio combining (H-S/MRC). If the signal fluctuation is measured at each branch output, the 
EFF is synonymous with the AF. In contrast, when the signal fluctuation is measured at the 
diversity combiner output, the EFF is equal to the square of the NSD of the instantaneous SNR 
at the combiner output. The term 'AF' was also confusingly used for diversity systems in some 
literature with a view to bridging the philosophy between characterizing physical channel fading 
and quantifying the degree of diversity effectiveness [42], [61], [62]. 

By definition, the efficiency of OSTBCs on reducing the severity of fading can be assessed 

by 

/JT^^ A /Var{7sTBc}\ 

EFFsTBc (dB) = lOlogio <^ -— -3 > 

UE{7stbc}) J 

= 101ogio{/t(||^||p)-l} (23) 
where k {\\H\\p) is the kurtosis of \\H\\p defined by 

E{[\\H\\^-E{\\H\\^}]'} 



n(\\H\ 



A 



" {E{[\\H\\^-E{\\H\\^}f})' 



2 ■ 



(24) 



In the second equality follows from the fact that the kurtosis is invariant with respect to 
translations of a random variable. Note that the minimum EFF is equal to —00 dB if there is no 
random fluctuation in the received signal. Also, the EFF is equal to dB for Rayleigh fading 
without diversity and hence, EFFstbc > dB means that the variation of the instantaneous SNR 
in each SISO subchannel is more severe than that in Rayleigh fading. 

1) Note on the Kurtosis of \\H\\p: The kurtosis measures the peakedness or flatness of a 
distribution [63]. It has been revealed that this normalized form of the fourth statistic of fading 
distributions plays a key role in the low-SNR behavior of the spectral efficiency in fading 
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channels [53], [64]. To proceed with deriving k (H/fHp) for double- scattering (tt-t, ^s, nR)-MIMO 
channels, we first define the following scalar quantity related to a correlation matrix. 

Definition 2 (Correlation Figure): For an arbitrary nxn correlation matrix the correlation 
figure of $ is defined by 

C(^>)^^ = ^tr($^) (25) 
tr (l^j ' 

where 1„ denotes the nxn all-one matrix. 

Note that ^ < C (^) ^ 1' where the lower and upper bounds correspond to uncorrected and 
fully correlated cases, respectively. The following Schur monotonicity properties hold for the 
correlation figure (the proofs are given in Appendix IIII-B I) . 

Property 1: Let $ be an ri x n correlation matrix. Then, the correlation figure ($) as a 
functional of the eigenvalues of $ is MIS, that is, if $ ^ then 

C($)<C(^). (26) 
Property 2: Let i = 1, 2, . . . , m, be x rij correlation matrices. Then, the product of 
correlation figures, Hi^i C i^i)^ a functional of the eigenvalues of <S)iLi^i^ is MIS, that is, 
if 



m m 



(27) 

i=l i=l 



then 

m m 



1=1 i=l 



Property 3: Let i = 1,2, ... ,m, be ra, x correlation matrices. Then, the sum of 
correlation figures, J2iLiC i^i)^ a functional of the eigenvalues of MIS, that 

is, if 

m m 
i=l i=l 

'"Similar to <25t . the correlation number was defined as i tr ($^) [54]. Wiiile tlie correlation figure and number are the 
second-order statistics of the spectra of a correlation matrix, normalized by those of fully correlated and uncorrelated matrices, 
respectively, the correlation figure is bounded by < ^ ($) < 1 for any correlation structure, as n ^ cxd. 
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then 

m m 

Ec(*^)<Ec(^0- (30) 

2=1 i=l 

The next theorem shows that /« (lliyjlp) depends exclusively on the spectra of spatial correlation 
matrices and is quantified solely by their correlation figures. 

Theorem 2: For double-scattering (nx, ns, nR)-MIMO channels, the kurtosis of \\H\\-p is 

k{\\H\\^) = C ($t) C (*r) + C (^-t) C (^s) + C ($r) C (*s) + 1. (31) 

Proof: See Appendix IIII-CI □ 

Example 2 (Spatially Uncorrelated Double Scattering): In the absence of spatial fading cor- 
relation (Tuc), we have 

K (||H||p) = + + + 1. (32) 

utUk utUs n^ns 

As compared with the i.i.d. case, the keyhole increases the kurtosis of the fading distribution 
in SISO subchannels by twice the reciprocal of the harmonic mean between the numbers of 
transmit and receive antennas, that is, — + — . 

Next, we show the Schur monotonicity property of k (||H||p). 

Corollary 1: Let 

j(T)A5l^e^^®^^^ (33) 

UtU^ n^ns ngnR 

for a spatial correlation environment T = ($t,^S5^r)- Then, the kurtosis of H-ffHp, as a 
functional of the eigenvalues of (T), is a MIS (or isotone) function, that is, if J (Ti) ^ (T2), 
then 

«:(||if||p;Ti)</s:(||H||p;T2). (34) 

Proof: It follows immediately from Theorem |2l and Properties |2l and |3] stating the fact that 
the product and sum of correlation figures preserve the monotonicity property. □ 

Corollary [T] implies that the less spatially correlated fading results in the less peaky fading 
distribution of each SISO subchannel. 
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2) Note on the EFF o/'Jstbc- From Theorem |2l and (l23b . it is straightforward to see that the 
EFFsTBc in double-scattering (nx, ^s, nR)-MIMO channels is given by 

EFFsTBc (dB) = lOlogio {C (^t) C (*r) + C (^t) C (*s) + C (^r) C (^s)} (35) 

from which we can make the following observations on the EFFstbc- 

• The EFFstbc as a functional of the eigenvalues of J" (T) is MIS, that is, 

EFFstbc (Ti) < EFFstbc (T2) (36) 

whenever (Ti) ^ (T2). This reveals that the less spatially correlated fading results in 
the less severe random fluctuations in equivalent SISO subchannels induced by OSTBCs. 

• In the absence of double scattering, ($s) is equal to zero and thus, the double scat- 
tering together with spatial correlation causes the EFFstbc to increase by the amount of 
C ($t) C (^s)+C (^r) C (^s)- In particular, the maximum increase in the EFFstbc is a sum of 
correlation figures of the transmit and receive correlation matrices, that is, ( ($t) + C (^r)' 
which eventuates when $s goes to be fully correlated or when the keyhole effect takes 
place. 

B. Low-SNR Capacity 

Recent information-theoretic studies show that the first-order analysis of the capacity versus 
the SNR fails to reveal the impact of the channel and that second-order analysis is required to 
assess the wideband or low-SNR performance of communication systems [53], [54]. In particular, 
it was demonstrated that the tradeoff between the capacity in bits/s/Hz and energy per bit required 
for reliable communication is the key measure of channel capacity in a low-SNR regime. In this 
regime, the capacity can be characterized by two parameters, namely, i) ^ . , the minimum 
bit energy per noise level required to reliably communicate at any positive data rate (where 
denotes the total transmitted energy per bit), and ii) Sq, the low-SNR slope (bits/s/Hz per 3 dB) 
of the capacity at the point ^ . . 

1) General Input Signaling: Before proceeding to study the low-SNR capacity achieved 
by OSTBCs, we first deal with the more general case of input signaling, assuming that the 
fading process is ergodic and coding is across many independent fading blocks without a delay 
constraint. 
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Theorem 3: Consider a general double-scattering {n-£,n^,n^)-M\MO channel 

Y = HX + W (37) 

where the channel matrix H is given by ^ at each coherence interval and the input signal 
X G C""^^^" is subject to the power constraint E{||X||p} = N^J^. Suppose that the receiver 
knows the realization of H, but the transmitter has no channel knowledge. Then, the minimum 
required ^ for reliable communication is 

^ loge2 

^Omin '^R 

and the low-SNR slope of the capacity is 



(38) 



2 

So = ^ ,^ X Tnr-^ ^ ,^ , ^ ,^ ^ ,^ ^ ,^ X bits/s/Hz per 3 dB. (39) 

C ($t) + C (*s) + C (*r) + C ($t) C (*s) C (^r) 

Proof: See Appendix IIII-Dl □ 
From Theorem |3l we can make the following observations. 

• The ^ . is inversely proportional to riR, whereas the double scattering and spatial fading 
correlation as well as the numbers of transmit antennas and effective scatterers do not affect 
this measure. Moreover, regardless of the number of antennas and propagation conditions, 
the minimum received bit energy per noise level required for reliable communication, 

"0 min 



is equal to 



^ =^R-jr =-1.59dB (40) 

' min ' min 



which is a fundamental feature of the channels where the additive noise is Gaussian [53, 
Theorem 1]. 

The low-SNR slope Sq as a functional of the eigenvalues of (T) is MDS, that is, if 
J{Ti) ^J{T2), then 

^o(Ti) >5o(T2) (41) 
where J" (T) is defined for the environment T = (^Tj^'sj^r) as follows: 

j(T)4il»giiffigaffi ^^«*^«*» . ,42) 

riT ns riR nxnsnR 
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Note that (BTTl follows from (l39b and Properties |2l and |3l This MDS property reveals that the 
low-SNR slope decreases with the amount of spatial correlation in contrast to the high-SNR 
capacity slope min (nT,n^,ns), which is invariant with respect to spatial correlation [21]. 

Example 3 (Dual-Antenna System): Consider n-r = = 2. In the presence of spatially 
uncorrelated double scattering, the low-SNR slope for general double- scattering (2, ns, 2)-MIM0 
channels is 

So = 2- + — bits/s/Hz per 3 dB (43) 



ns 4^ 

which is bounded by 8/9 < 5*0 < 2. The lowest and highest slopes are achieved when ns = 1 
(keyhole) and ns = oo (i.i.d.), respectively. 

2) OSTBC Input Signaling: We now turn attention to the low-SNR behavior of the capacity 
for double- scattering (nx, ns, nR)-MIMO channels employing OSTBCs. 

Theorem 4: Consider a double-scattering (tt-t, ns, nR)-MIMO channel 

Y = HQ' + W 

where the channel matrix H is given by ^ at each coherence interval and the OSTBC Q is 
subject to the power constraint E|||^||p} = N^V. Then, the OSTBC achieves the minimum 
required ^^^^ same as that without the orthogonal signaling constraint 

i^^b^™^ log, 2 



(44) 



^0 min '^R 

and the low-SNR slope of the capacity 

27? 

-^n™^ = s , , , , bits/s/Hz per 3 dB. (45) 

Proof: See Appendix UlTEl □ 

From Theorem |4l we can make the following observations in parallel to IIV-B.II 

• As compared with the general case, the use of OSTBCs does not increase the minimum 
required ^ for reliable communication in MIMO channels. 

• The low-SNR slope Sq^^'~ as a functional of the eigenvalues of (T) is MDS, that is, if 
J'(Ti) ^J{T2), then 

5o'™'^(Ti)>5o^™^(T2). (46) 
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In contrast, we see from (11591 ) that the high-SNR slope of the capacity is equal to TZ, which 
does not depend on spatial correlation and double scattering. 

Example 4 (Alamouti's Code): Consider = n-g^ = 2. In the presence of spatially uncor- 
related double scattering, the low-SNR slope for Alamouti's code with two receive antennas 
is 



-<STBC " ' - 14 



^^ibc = r . ( 1 + ^ . ^ ) bits/s/Hz per 3 dB (47) 

5 V ns 5J 

which is bounded by 8/9 < Sl™^ < 8/5. 

In Fig. |7J the capacity (bits/s/Hz) versus . and its low-SNR approximation are depicted 
with and without the signaling constraint of the OSTBC Q4 in double-scattering (4, 20, 4)-MIM0 
channels with exponential correlation $x = = ^4'^'' (0.5) and $3 = $20* (0-5). For the 
OSTBC Ga, the low-SNR approximation is remarkably accurate for a fairly wide range of ]Vomin' 
whereas there exists some discrepancy between the Monte Carlo simulation and the first-order 
approximation for the general input signaling — approximately 11% difference at . = dB, 
for example. In this scenario, the low-SNR slopes are 1.26 and 2.46 bits/s/Hz per 3 dB with 
and without the OSTBC input signaling constraint, respectively. Thus, the use of the OSTBC 
^4 incurs about 49% reduction in the slope. This slope reduction is much smaller than that in 
a high-SNR regime: the high-SNR slope for the OSTBC ^4 is 7^ = 3/4 and the corresponding 
slope for the general signaling is equal to min (nx, n^,ns) =4 bits/s/Hz per 3 dB [21]. 

V. Conclusions 

We investigated the combined effect of rank deficiency and spatial fading correlation on the 
diversity performance of MEMO systems. In particular, we considered double-scattering MIMO 
channels employing OSTBCs which use up all antennas to realize full diversity advantage. 
We characterized the effects of double scattering on the severity of fading and the low-SNR 
capacity by quantifying the EFF and the capacity slope in terms of the correlation figures of 
spatial correlation matrices. The Schur monotonicity properties were shown for these performance 
measures as functionals of the eigenvalues of correlation matrices. We also determined the 
required scattering richness of the channel to achieve the full diversity order of nxriR. Finally, 
we derived the exact SEP expressions for some classes of double scattering, which consolidate 
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the effects of rank efficiency and spatial correlation on the SEP performance. On account of the 
generality of channel modeling, the results of the paper are substantial enough to encompass 
those for well-accepted existing models (e.g., i.i.d./spatially correlated/keyhole MIMO channels) 
as special cases of our solutions. 

Appendix I 

Majorization, Schur Monotonicity, and Correlation Matrices 

We use the concept of majorization [65]-[69] as a mathematical tool to characterize different 
spatial correlation environments. Using the majorization theory, the analytical framework was 
established in [52] to assess the performance of multiple-antenna diversity systems with different 
power dispersion profiles. In particular, monotonicity theorems were proved for various perfor- 
mance measures such as the NSD of the output SNR, the ergodic capacity, the matched-filter 
bound, the inverse SEP, and the symbol error outage. The notion of majorization has also been 
used in [18], [36], [70] as a measure of correlation. In this appendix, we briefly discuss the basic 
properties of majorization and Schur monotonicity. 

A. Majorization and Correlation Matrices 

Given a real vector a — (ai, a2, ■ ■ ■ , e M'*, we rearrange its components in decreasing 
order as a[i] > op] > • • • > a[„]. 

Definition 3: For a = (oi, a2, ■ ■ ■ , a™)^, b = {bi,b2, ■ . . , bn f G M", we denote a -<b and say 
that a is weakly majorized (or submajorized) by h if 

k k 

A; = l,2,...,n. (48) 

i=l i=l 

If Yll=i = Si=i h holds in addition to a 6, then we say that a is majorized by 6 and denote 
as a ^ 6. 

For example, if each Oj > and Y^^=i — then 

^ (ai,a2,...,a„)^^ (n,0,...,0)^. (49) 

Of particular interest are the majorization relations among Hermitian matrices in terms of their 
eigenvalue vectors to compare different spatial correlation environments. A Hermitian matrix A 
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is said to be majorized by a Hermitian matrix B, simply denoted hy A ^ B, if X (A) ^ A (B) 
where A (■) denote the vector of eigenvalues of a Hermitian matrix. For example, the well-known 
Schur's theorem [68, eq. (5.5.8)] on the relationship between the eigenvalues and diagonal entries 
of Hermitian matrices can be written as 

AoI^^A for Hermitian A e C'^^'^ (50) 

where o denotes a Hadamard (i.e., entrywise) product. One of the most useful results on the 
eigenvalue majorization is the following theorem. 

Theorem 5 ( [67, Theorem 7.1]): A linear map C : C"^" ^nxn called positive if 
£ (A) > for j4 G C"^" > and unital if C {In) = In- It is said to be doubly stochastic 
if £ is a unital positive linear map with the trace-preserving property, i.e., iiC{A) = ii{A), 
MA E C"^". Let A e C"^" be Hermitian and £ be a doubly stochastic map. Then, 

£{A)^A. (51) 

Recall that the Schur product theorem [68, Theorem 5.2.1] says that the Hadamard product of 
two positive semidefinite matrices is positive semidefinite. Therefore if $ G C"^" is an arbitrary 
correlation matrix and define £ (A) = Ao^, then £ is obviously a doubly stochastic map on 

Corollary 2: Let A G C"^" be Hermitian and $ G C"^" be a correlation matrix. Then, 

Ao^^A. (52) 

In fact, this result was first given in [69, Corollary 2] without using the notion of doubly 
stochastic maps. From Corollary |2l we can obtain the eigenvalue majorization relations for 
the well-known correlation models — constant, exponential, and tridiagonal correlation — which 
have been widely used for many communication problems of multiple- antenna systems (see, 
e.g., [21]-[23], [49], [54], [71]). 

Example 5 (Constant, Exponential, and Tridiagonal Matrices): The nth-order constant, expo- 
nential, and tridiagonal matrices with a coefficient p, denoted by (p), ^[^^ (p), and (p) 
respectively, are n x n symmetric Toeplitz matrices of the following structures: 
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(P) 



^^r^ (p) 



^i'^ (p) 



-L p 


p ■ ■ 


■ P 


p 1 


p ■• 


■ P 


.p p 


p ■ ■ 
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" 1 
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1 
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pn-2 




ri P 






p 1 


P 




p 


1 


P 



(53) 



P 



("-1)1 

.(«-2) 



(54) 



1 P 
P 1. 



(55) 



Note that (p), (p) with p G [0, 1] and ^^^h p G [O, 0.5/ cos are correlation 

matrices, since they are positive semidefinite for such values of p. Let < pi < p2. Then, since 

(Pl)=<'^^ (P2)o<'=) (f^^ 
^i'^ (Pl)=<'^ (P2)o<'=) 

it follows from Corollary |51 that 



^i^^ (Pl) ^ ^i^^ (P2) 
^i'^^ (Pl) ^ ^i'^^ (P2) 
(Pl) ^ (P2) • 



(56) 
(57) 
(58) 



Remark: If < pi < p2, then and are positive semidefinite. Hence, the 

majorization relations (l56t - (l58t hold, although each matrix itself is only Hermitian but may not 
be positive semidefinite. 

B. Schur Monotonicity 

The concept of majorization is closely related to a MIS (or MDS) function. If a function 
/ : (a subset of) M" -> M satisfies / ( ) < f {bi, ■ ■ ■ ,bn) whenever a ^ b, then / is 
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called a MIS (or isotone) function on (a subset of) W\ The following theorem gives a necessary 
and sufficient condition for / to be MIS. 

Theorem 6 (Schur 1923): Let I C M and / : I" ^ R be continuously differentiable. Then, 
the function / is MIS on I" if and only if 

/ is symmetric on I"' (59) 

and for all i j, 

df df 



> Va e r. (60) 



Note that Schur' s condition (l60b can be replaced by 

df df 



> Va G r (61) 



dai da2 

because of the symmetry. If / is MIS on I", then — / is a MDS function on I". 

Appendix II 

Some Statistics Derived From Complex Gaussian Matrices 
This appendix gives useful results on some statistics derived from complex Gaussian matrices. 

A. Preliminary Results 

Lemma 1: Let ~ -Afm.n (Omxn, ^jt), k = 1,2, . . . ,p, he statistically independent com- 
plex Gaussian matrices and 

X = [Xi X> ■ ■■ Xp\ ~ J\fm,np (Omxnp, S, 0fc=i ^ k) ■ (62) 

Then, for A e C™^'" > and 5 = 0^^^ Bk, Bk G C^""" > 0, we have 

p 

E {etr {-AXBX^) } =Y[det (/„„ + ® ^kBk)~' ■ (63) 

k=l 

Proof: Since AXBX^ = J2l=i AX kB kXl, we have 

p 

E{etr {-AXBX^)} = JJe^, {etr (-ylX.BfcXy } . (64) 

k=l 
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Therefore'^ 

Ex, [etT{~AXkBkXl) 

= cj eii{-AXkB,X\-i:-^X,^,'Xl)dXu 

= CkJ^ exp [- (vec(Xi))^ { (E^ ® ^,)"' + ® 5,} vec{Xl) 

= Cfc7r'""det|(E^®^fc)~'+A^®Bfc}"' 
= det (7™„ + AE®^fc5fc)~^ 

where = 7r^™"det (E)^"det (^fc)^™. Combining (l64b and (l65l) complete the proof. 
Lemma 2; Let X ~ AT^^n (0„xn,E,^). Then, for A,B e C"''", we have 

E {etr (XU + B^X) } = etr (EA^B^) . 
Proof: Let Mi and M2 be m x n matrices such that 
tr(XU + B'X - E-^X^-^X^) 

= tr(E 'M,^ 'M^) + tr |-E ' (X - Af,) (X - M-.f^ . 

Then, since 

/ etr{ -E-' (X - Ml) (X - M-. fjdX = tt'"" det (E)" det , 

we get 

E{etr(xU + E"^X)} = etr(E"^Mi^"^M5) . 
By comparing both the sides of (l67b . we have 

Ml = EA^ 
M2 = E5^. 

Finally, substituting dTDI) and dTTT) into ( 15^ completes the proof. 

"if X = (Xij) is an m X n matrix of functionally independent complex variables, then 

m n 

dX = Y{Y{ '^^^^J dQX^■j. 
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Lemma 3: Let X ~ N'm,n (M.E.^). Then, the characteristic function of X is 



<^x {Z)^E{exp[jmT (XZt)]} 



exp 



(72) 



where j = y/—T and Z G C"^" is an arbitrary matrix. 
Proof: Let Xi ~ AL,„ (0„xn, E, ^). Then, 

<l>x{Z) =exp [j3fttr(MZ^)] ■ E {exp [j3ft tr (XiZ^)] } . (73) 

Since 

3ft tr (XiZ^) = ^ tr [Z^Xi + XjZ) , (74) 
it follows from Lemma El that 

E {exp [j3fttr (XiZ^)] } = etr (^-^EZ^Z^^ . (75) 

Combining (l73t and (l75t completes the proof. □ 
We remark that Lemma |3] is a counterpart result of the real case in [72, Theorem 2.3.2]. 

B. Hypergeometric Functions of Matrix Arguments 

The hypergeometric functions of matrix arguments often appear in deriving the distributions 
and statistics of random matrices [72]-[76]. In parallel to the hypergeometric functions of a 
scalar argument, the hypergeometric functions of one or two matrix arguments can be expressed 
as an infinite series of zonal polynomials:^^ 

A:=0 K I- ^ "^Jk 

'^Zonal polynomials of a symmetric matrix were introduced in [73] using group representation theory. In parallel to a real 
matrix argument, zonal polynomials of a Hermitian matrix were defined in [74] as natural extension of the real case. Those 
polynomials are homogeneous symmetric functions in the eigenvalues of matrix argument and can be constructed in terms 
of homogeneous symmetric polynomials such as monomial symmetric functions, elementary symmetric functions, and Schur 
functions [77]. 
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F(n)/, , - J,. A ON ["iJ. • • • {A) q {B 



A:=0 



(77) 



with Hermitian A G C"^" and 5 G C"^". In (EH) and dTTl), k = [h, k2, . . . , K) denotes a 
partition of the nonnegative integer k such that ki > k2 > ■ ■ ■ > kn > and Yl^=i = ^' 
[a]^ is the complex multivariate hypergeometric coefficient of the partition k, [74, eq. (84)], and 
Ck{-) is the zonal polynomial of a Hermitian matrix [74, eq. (85)]. Although these functions are 
of great interest from an analytical point of view, the practical difficulty lies in their numerical 
aspects. The determinantal representation for the hypergeometric function of two Hermitian 
matrices [76, Lemma 3] settles this computational problem and has been widely used in the 
literature of multiple- antenna communication theory (see, e.g., [22], [23], [55], [56]). However, 
[76, Lemma 3] is valid only for the case of two matrix arguments with the same dimension 
and the distinct eigenvalues. In the following lemma, we generalize [76, Lemma 3] for the case 
that two matrix arguments have the different matrix dimension and the eigenvalues of arbitrary 
multiplicity. 

Lemma 4 (Generic Determinantal Formula): Let A G C™^™ and S G C"^", m < n, he 
Hermitian matrices with the ordered eigenvalues Ai > A2 > • . • > and ci > a2 > . . . > a„, 
respectively. Given Oj, bj G C where i = 1,2, . . . ,p and j = 1,2, . . . , q, define 

Kpg (x) = pFg (ai - n + u, . . . , ap - n + iy;bi - n + u, . . . ,bg - n + u; x) (78) 

''"^ nLi(«.-^+i). ^ ^ 

where v is an arbitrary nonnegative integer, (a)^ = a (a + 1) ■ ■ ■ (a + n — 1), (a)^ = 1 is the 
Pochhammer symbol, and pFq (ai, 02, . . . , a^; 61, 62, • • • , 2) is the generalized hypergeometric 
function of scalar argument [78, eq. (9.14.1)]. Then, 



pFf ) (ai, . . . , ap; 61, . . . , 6,; A, S) ■ J] (A,- - A,) 

i<j 



p,q 



3^2 



det (A)" det {[Z(^n),i Z(n),2 ■■■ 



(80) 
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with 

m 
i=l 

where yu = (3^fc,ii) and 2(^1^^^ = {Z(^i)^k,ij), I < n, k = 1, 2, . . . , ^ (S), are m x (E) and 
I X r/c (S) matrices, whose (i, j)th entries are given respectively by 

yk,ij = T • Hp'fg (^if^(fe>) (82) 
2(l•,,k,^J = {i-J + (^\ky (83) 

In particular, for q-Fq"^ (A, 2), -ft'"^" in dSB and the {i,j)th entry of yk in (US) reduce to 

m 



^oT = n(^-^)' (84) 
3^Mi = Ar'e^'^^w. (85) 



Proof: Let us dilate the m x m matrix A to the n x n matrix A © On-m by affixing 
zero elements. Then, this augmented matrix A © On-m has the eigenvalues Ai, A2, . . . , and 
[n — m) additional zero eigenvalues. Note that zonal polynomials depend on its Hermitian matrix 
arguments through Schur functions in the eigenvalues of matrix arguments [74]-[77]. Since Schur 
functions are invariant to augmenting zero elements [79], it is easy to show that 

a(A©0„_™) =a(A). (86) 

Let Am+i, Am+2, . . . , A„ be [n — m) additional zero eigenvalues and denote the left-hand side of 
(ISUt by LHSj8o| for convenience. Then, it follows from (15^ and [76, Lemma 3] that 

^ det^ J^pi (A.^,)) m 
LHS|8o| = i^p/g" prn " ryT- -r ■ TT (Aj - A^) . (87) 

From a computational point of view, (1571) presents numerical difficulty since the Vandermonde 
determinant YYi<j (Aj ^ Ai) or YYi<j i'^j ~ '^i) becomes zero when some of the Aj's or o-j's are 
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equal. This can be alleviated by using Cauchy's mean value theorem (or L'Hospital's rule): 

det {n;:i{Ka^)) m 

LHSiHol = K;," lim lim t^^^ t^t r ■ TT (^J' ' (^8) 



where a ^ a means that 



ri(S) 

CT/ 



i \ri(S)+r2{E) 



{^»}i=„_r,(5:)(E)+l ^ %(S)>- 

Let n-dimensional vectors u {z) and v (z) be 

n (z) = {n;:l {a,z) , n;:l (a^z) , . . . , n;:l (89) 

v{z) = {l,z,...,z^-') (90) 

and let u^'^^ (z) and v^''^ (z) be the kth derivatives of u (z) and v (z) with respect to z, respectively. 
Note that the jth components Uj'^ (z) and f j'^^ (z), j = 1, 2, . . . , n, of u^^^ {z) and t;'^''') (2;) are 
given respectively by 

k 

^ i^) = %- Ki^' (^.-) (91) 

vf\z) = {3-k),z^-^-' (92) 

where follows from the differentiation identity of [80, eq. (7.2.3.47)]. Then, taking the limits 

on Afc's, we get 

lim - = (93) 

{Aa^.™,.^o nr<,(A,-A.) det([^^]) ^ ^ 

with the (n — m) X n matrices 

July 5, 2006 DRAFT 



30 



REVISED FOR PUBLICATION IN THE IEEE TRANSACTIONS ON INFORMATION THEORY 



Ua = (f/A,.,) 



u 



(0) 
(0) 

(n— m— 1) 

(0) 
(0) 



(94) 



(0) 



(95) 



V 



(n—m—l) 



(0) 



and the m X n matrices Ub = {'H'p'l (AjCTj)) and Vb = (A] ). From (|9T1) and (l92b . it is easy 
to see that the («, j)th entries of Ua and Va are given respectively by 



v^-'^ (0) 



Va., = (0) 



n,i—l 

Xp,q 



(2-1)!, if 2 = J 
, otherwise. 



Now, using the result on the determinant of a partitioned matrix 

^ = det (A) det {D - CA'^B) , if A is invertible, 



det 





A 


B 


( 


C 


D 



we have 



det 



Va 
Vb 



1)! ■ det 



1=1 



/ 


-^n—m 




-m+1 


■ A^i- 


\ 




\ n—m 
^2 


\ n- 
^2 


-m+1 


\ n—l 
■ ^2 




\ 


\ n—m 
L^m 




-m+1 


\n—l 
■ J 


/ 



Z=l A;=l 

Hence, combining (l88t . (l93t . and (l99t gives 



LHS 



T/'m.n 



det 



lim 



det (A)"-" nr<, - 



(96) 
(97) 



(98) 



(99) 



(100) 
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where C/a = (o"] ^) is the (n — m) x n submatrix of the Vandermonde matrix of ai, (J2, . . . , cr„. 
Using similar steps leading to (l93l) . we obtain 



det 



lim 



det 



•^•{n-m),! ^{n-m),2 



yi 



y2 



^{n-m.),Q(£) 



n 



n ( 

i<3 



det ( \Z(ri 



),1 ^(n),2 



]) 



(101) 



where the (i,j)th entries of m x (E) matrices 3^^ and / x (E) matrices fe, / < n, 
= 1, 2, . . . , (S), are given by (l82b and (l83t . respectively. Finally, substituting (IIOII) into 
(1 1001) completes the proof of the lemma. □ 

As a by-product of Lemma IH we obtain the following determinantal formula for the hyper- 
geometric function of one matrix argument. 

Corollary 3: If E = /„ in Lemma |4l then we have 



{a„ . . . , ap; 61, . . . , 6,; A) ■ n (A,- - A,) = det (A^^7^^;-'"+^■ (A,)) . 



(102) 



Proof: The result follows immediately from d^St and Lemma|?|with g (E) = 1, ti (E) 



and (T^-^^ = 1 



n, 



□ 



C. Some Statistics 

Lemma 5: Let X ~ Afm,n (0^xn,S,'^)- Then, for A G C™^" > and 5 G C"^" > 0, the 
/cth-order cumulant of ti [AXBX'') is 

Cum,{tr (AXBXt)} 4 ^ In 0,,(^;,^^t) (s) 

= (A;- 1)! tr{(AE)*^}tr{(^B)'^} (103) 

where 0tr(AXBXt) («) = E{etr (-sAXBX"^) } is the MGF of tY(AXBX^). 
Proof: Since 

tT{AXBX^) = {vec{X^)y [A^ ® B) vec{X^) 
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is a quadratic form in complex Gaussian variables, whose characteristic function has been 
reported in [81], it can be readily shown that 



r(AXBX^) (s) = det + s(S^ ® ^) [A^ ® B) } 
= det (Imn + sAT. (g) ^S)"^ . 



-1 



Therefore, 

^ln0tr(Axsxt)(s) 



-1)^'(A; - l)!tr 



sAE ® ^B)-^ (AE ® 



(104) 

(105) 

□ 



Hence, we obtain the result ( 11031 ) from (11051) with s = 0. 

We remark that the cumulants, except for the first-order cumulant, are invariant with respect 
to translations of a random variable. The first and second order cumulants are the mean and 
variance of the underlying random variable, respectively, and other higher-order statistics can 
also be obtained from general relationships between the cumulants and moments. Lemma |5l 
reveals that all cumulants of tii^AXBX'') as functional of the eigenvalues of AT, and '^B are 
MIS. 

Lemma 6: Let X ~ Afrn,n (0„xn,S,^). Then, for A G C™^'" > and 5 e C"^" > 0, we 
have 



EUr 



[AXBX^Y]] =tr^ (AE)tr{(*B)'} +tr2 ($5) tr{(AE)'}. 



(106) 

Proof: We first start with the characteristic function of S = {Sij) = A^^'^XB^^'^. Let 
E = [tij] = A^I'^Y.A^I'^ and ^ = {^ij) = B^^^^B^^^. Then, 



^siZ) =E|exp [j3?tr(A'/'X5'/'Zt) 

= $x (a'/'zb'/' 



(a) 



etr 



4 



(107) 



where (a) follows from Lemma |3] and 



m m n n 



(108) 



i=l p=l q=l j = l 
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It follows from the characteristic function $5 (Z) in (11071) that 



E { Si^j^ S'*,j2 5'j3j3 S*^j^ } 



1 



d<^s {Z) 



8^3 jZ) ^ ^^3 jZ) 
dnZi 



UJ4 



d'^Z 



14]4 



with 



(Z) 



^2 {Z) 



^3 {Z) 



d^Z, 



dyi {Z) 

[Z) 

d^Zi.^j^ 



J 



d^i {Z 



d^Z„ 



d^Zi^j^ 
dip2 {Z) 



d'^Z 



«3J3 



Using (I109I) . we obtain 



Ex tr 



m n n m 



XI X] X] ^ {^iP^jp^jg^iq} 

i=l p=l q=l j = l 



m n n m 



^Sjj^ppSjj^gg + SM^p<jSjj\i' 

1=1 p=l g=l j=l 



(109) 



(110) 

(111) 
(112) 



(113) 



□ 



from which (11061) follows readily. 

Theorem 7: Let Xi ~ A/'m^p (Omxp, Si, ^1) and X2 ~ ATp .„ (Opxn, ^2, ^2) be statistically 
independent complex Gaussian matrices. Then, 

Ex,,x. [iT^{X,X^X\x\)] 

= tr (E?) tr^ (^iE2) tr (^2) + tr (E?) tr^ (^2) trK^iEs)'} 

+ tr^ (El) tr{(^iE2)'} tr {^l) + tr^ (Ei) tr^ (^iE2) tr^ (^2) (114) 
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and 



E 



tr 



= tr^ (El) tr^ (^iE2) tr {^l) + tr^ (Ei) tr^ (^2) tr{ (^lEs)'} 
+ trK^iSs)'} tr {^D tr (E?) + tr^ (^^Es) tr^ (^2) tr (E?) 

Proof: Using the first two cumulants from Lemma |5l we get 

Ex„x. [ir'{X,X.,X\X\)] 

Ix, {tr (E?) tr [{X.,X\^, f] + tr^ (Ei) tr^ {X^X\^,) } 



wliere it follows from Lemma |6l that 



Ex. <i tr 



{X,X\^,f 
and from Lemma |5l that 

Ex, {tr^ {X^X\^,) } = tr { (^iE2)' } tr {^l) + tr^ (W1E2) tr^ 

Combining (II 16l) - (ll 181) yields the desired result (11141) . 
Similar to (11161) . we have 



Exi,x. <! tr 



(-^1-^2-^2-^1) 

Ex, {tr2 (Ei)tr [(X.^t^,)'] + tr^ (X2X5^i) tr (E?) } 



(115) 



(116) 



tr2(^iE2)tr(^2) +tr2(^2)tr{(^iE2)'} (117) 



(118) 



(119) 



□ 



From (II 17l )- (ll 191 ). we obtain the desired result dl 151 ). 

Theorem 8: Let X ~ N'm,n (Omxn; -^n)j m < n, and (Ji, o"2; • • • ? be the eigenvalues of 
E in any order. Then, the joint pdf of the ordered eigenvalues Ai > A2 > . . . > > of a 
central complex Wishart matrix XX^ ~ VVm {n, E) is given by 



Pa(Ai,A2,...,A„) =^-Met([Gi G2 



l<«,7<m ■' ' -M-J. 



(120) 



where 



DRAFT 



July 5, 2006 



SHIN AND WIN: MIMO DIVERSITY IN THE PRESENCE OF DOUBLE SCATTERING 



35 



(121) 



and Gk = {Gk,ij) and = [Bk,ij), k = 1,2, . . . , g (E), are m x (E) matrices, whose («, j)th 
entries are given respectively by 



Gk,ij 



■ I 1 \ n—i+j 



(122) 
(123) 



Proof: The joint eigenvalue density px (Ai, A2, • • • , Am) is given by [74, eq. (95)] in terms 
of the hypergeometric function of matrix arguments. To render this joint pdf more amenable 
to further analysis and computationally tractable, we apply Lemma lU to [74, eq. (95)], which 
results in (11201) after some algebra. □ 

Note that (11201) is valid for any covariance matrix E with the eigenvalues of arbitrary multiplic- 
ity and hence, generalizes the previous determinantal representation for the joint eigenvalue pdf of 
Wishart matrices. If E = in Theorem HI all of the eigenvalues are identically equal to one and 
hence, with g (E) = 1, ri (E) = m, and cr^^^ = 1, (11201) reduces to [22, eq. (6)]. Furthermore, if all 
the eigenvalues of E are distinct, then, with g (E) = m and ti (E) = T2 (E) = . . . = (E) = 1, 
(fT20b reduces to [22, eq. (18)]. 



Theorem 9: Let X ~ N'rn,n {Qmxn, Im,'^), m < n, A E C"^" be Hermitian positive definite, 
and Pi, P2, ■ ■ ■ , Pn be the eigenvalues of A^^'^^A^^'^ in any order. Then, the joint pdf of the 
ordered eigenvalues Ai > A2 > • • • > Am > of a matrix quadratic form XAX^ is given by 



Px (Ai, A2, . . . , A^ 



det ^ 


V(n-m),l V(„_m),2 

Qi Q2 


■ ■ ■ "^(n-m),e(Ai/2*Ai/2) 


) 


i^oT^et (A^r det ( 


V(„),i 


V(n),2 ■ ■ ■ ^(n),g{A'/^^A'/^) 


) 



det (A;.-i) (124) 

l<i,j<m ^ ■' ' 



where = {Qk,^^) and V(«),fe = (V(/),fc,,j), I < n, k = 1,2, . . . , g{A^/^^A^^^), are 
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Tfc (A^^^^^^^^) and / X Tk[A^^'^^A^^'^^ matrices, whose (i, j)th entries are given respectively by 



Qk,j = Ar^e-"'/''« (125) 

j-lP{k) ■ 



V(o,.,, = (-1)^'^' i^-J + P7^'. (126) 



Proof: Let S = XAX\ then S ~ Qm „(A,/m,^) is a positive-definite quadratic form 
in the complex Gaussian matrix [21, Definition 11.3]. Using the pdf [23, (2)], we can write the 
joint eigenvalue pdf of S in the form 

_m(m— 1) r 

(Ai, A2, . . . , A„) = - / PS (UDW) TT (A, - A,)' [dU] 

r„ (n) r„ (m) ^^^^ f^-^ 

(127) 



where D = diag (Ai, Aa, . . . , A„), (a) = n'^i^-^)/^ T (a - i) with 3? (a) > m - 1 is 

the complex multivariate gamma function and T (■) is the gamma function. In (I127I) . U (m) = 
[U : UU'^ = Im] is the unitary group of order m and [dU] is the unitary invariant Haar measure 
on the unitary group U (m) normalized to make the total volume unity. Similar to Theorem [H 
we obtain the desired result (11241) applying Lemma |4| to (11271) . □ 

Definition 4 ( Characteristic Coefficient): Let A be an n x n Hermitian matrix with the eigen- 
values ai,a2, ■ ■ ■ ,an in any order. Then, the {i,j)th characteristic coefficient Xij{A), i = 
1,2, . . . , g (A), j = 1, 2, . . . , Tj (A), is defined as a partial fraction expansion coefficient of 
det (/„ + ^Ay^ such that 

det (7„ + eA)-^ = n(l + e«(.))"''^''^ 

1=1 

s(A)n{A) 

= EE'^'.:'(^)(l + ^"»)"' (128) 
i=i j=i 

where ^ is a scalar constant such that In+^A is nonsingular. The {i, j)th characteristic coefficient 
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Xi j (A) can be determined by 



{A) 



1 



a 



E 



(«) A;i+fc2+...+A:p(^)=roij' /=! 
fc,e{0,N} for Vi^i 
fci=0 



a 



(0 



(129) 



where = (A) — j. 

Note that the characteristic coefficients are invariant with respect to the constant and only 
a function of the spectra of A. In addition, it can be seen from (11281) with = that the sum 
of all the characteristic coefficients is equal to one. By definition, we have 



'^Ij {In J 



0, j = l,2,...,n-l 

1, j = n. 



(130) 



Example 6 (Constant Correlation Matrix): Consider a constant correlation matrix ^^^\p). 
Since the eigenvalues of $1^^ (p) are 1 + (n — 1) p and 1 — p with n — 1 multiplicity, it is 
easy to show that the characteristic coefficients of (p), p G (0, 1), are 



'^iA^^Hp)) = ( , \ (131) 
\l - p + np J 

X, , (p)) = . ( ^ (132) 

where j = 1, 2, . . . , n — 1. 

Theorem 10: Let X ~ ATm^^ (Omxn; ^j-'^n)? rn < n, and (Ji, o"2; • • • ? be the eigenvalues of 
E. Let Abe a u X u positive-semidefinite matrix with the eigenvalues ai, a2, . . . , Then, for 
^ > 0, we have 



E 



|det (7„, + eA®XXt) '} =^-Met ([Oi ■■■ 0,(e)]) (133) 
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where A is given in (11211) and = i^k,ij), k = 1,2, . . . , g (S), are m x (E) matrices whose 
(i, j)th entry is given by 



(,{A)r^{A) 

^k,ij 

p=l q=l 



5^ 5^ { A',,, (A) 4--+^+^-^ (n - m + ^ + J - 2)! 

X 2Fo{n-m + i+j - l,q; -^a{p)Cr{k))] (134) 



where g (A) is the (p, g)th characteristic coefficient of A. 
Proof: From Theorem IH we have 



E 



p=i 



J-N -Tp(A) 



■■■ n + ea(p>Afc)""''%A(Ai, A2, . . . , Xm)dXidX2 ■■■dK 

^ k=l p=l 



0<A™,<...<Ai<oo 



X det (\Gi G2 ■■■ det (A;-^) dX^dX^ ■■■dX^ 

^^A-^det{[ni ^2 ••• (135) 

where (a) follows from the fact that the integrand is symmetric in Ai, A2, . . . , A^ and (6) follows 
from the generalized Cauchy-Binet formula [22, Appendix], [23, Lemma 2], yielding the (i, j)th 
entry of m x (S) matrices Qk, k = 1,2, . . . , g (S), as 



^k,ij 



poo (>(A.) 

/ r[(l + e"(p)A)"''"^'^^A"-"^+'+^-2e"^/"«dA. (136) 
Jo 1 



'0 p=l 

Using a partial fraction decomposition, (11361) can be written as 



'IJV— ^ poo 

^k,j=J2J2^P''^^^'^ / (l + e«(p)A)"V-"^+^+^-V'^^)rfA (137) 
p=i 13=1 
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where the characteristic coefficients Xp g (A) is given by (I129I ). We complete the proof of the 
theorem by evaluating the integral in (11371) with the help of the following integral identity: 



(138) 



poo 

/ (1 + axf-^ x"-^e-"/^dx = 6" (n - 1)! (n, -/x + 1 ; -ab) 
Jo 

where a,b > 0, n E N, and fi E C. □ 
Corollary 4: Let X ~ Afm,n (Omxn,S,/„), m < n. Then, for G N, we have 

E{det(/„ + eXXt)-n = /"^^.^^ (139) 

where il = (i^ij) is the m x m Hankel matrix whose («, j)th entry is given by 

n^j = (n-m + i+j - 2)! 2F0 (n - m + i + j - 1, u ; . (140) 

Proof: It follows immediately from Theorem [10| with E = 7^, A = I^, g (E) = 1, 

Ti (E) = m, a^-^^ = 1, ^ (A) = 1, Ti (A) = u, and a^-^^ = 1. □ 

Theorem 11: Let X ~ Afm.,n (Omxn, cTi, i = 1,2, ... ,m, and j = 1, 2, . . . , n, be the 
eigenvalues of E and ^, respectively. Then, for ^ > 0, we have 

e(E) £.(*)rp(S)r,(«) 

E|det(7„ + eXXt)-'} = ^ ^ A',,(E)A',,, W 2Fo(z,j;-ea<,^^<,^) (141) 

p=l qr=l i=l j = l 

where Xp^i (E) and Xgj (^) are the {p,i)th and (g,j)th characteristic coefficients of E and ^, 
respectively. 

Proof: It follows from Lemmas [T] and El that 

det {Im + e^Xt)"' = E,^ jetr (-^^^^ii/lx) } 

= E,^,,, {etr {^y\Xy, - v\X^v,) } (142) 

where ~ A/'m,! (0^x1,-^™, 1) and t/g ~ -Mi,! (Onxij-'^n, !)• Denoting the left-hand side of (11411) 
by LHS(i4i] and using (11421) . we have 

LHSen = {Ex {etr (^^2/1^ - X^ViV^^ } } 

= ¥.y^.y^ jexp [-iv\i:y,y\^y^^ } . (143) 
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Now, introducing a delta function to decouple the expectations for and t/g iii (11431) yields 

f POD POO / + \ 1 

L J —oo J ~oo ) 

'■ ¥.y^ |etr(-j^E?/i?/l) | E^,^ ^eii {-^z^y^yl) | dudz 

/CO poo 
/ e^^" det {Irn + JcuE)-^ det (J„ + ^z^)'^ dujdz 
■oo J —oo 



(a) _1_ 

27r 



1 

(J ^ 
27r 



/•oo POO 



-oo J —oo 



e(S)e(*)rp(E)r,(*) 
p=l 5=1 i=l j=l 



oo POO 



X / e^-^(l+ja<p^^)"Xl + e^(5>^)"'^^^4 (144) 



oo J — oo 



where (a) is obtained by replacing the delta function with its Fourier representation, (6) follows 
from Lemma [H and (c) is obtained from Definition |4| Using the integral identity, for a > 0, 

£ G N, and 2 G M, 

e^-' {l+jauj)-'du = ""^^ (1 + sign {z)) , (145) 

(11441) can be written as 

LHSoai = E E E E ^.'^^i r' / (1 + »,„.)-'.-e-/'«<i.. (146) 

p=i i=i j=i ^{p) y^- Jo 

Finally, we obtain the desired result (11411) by evaluating the integral in (11461) with the help of 
(fT38l). □ 



Appendix III 
Proofs 

A. Proof of Theorem \T\ 

We first prove Theorem [T] for Af-ary phase shift keying (M-PSK) signaling. The SEP of the 
OSTBC with M-PSK constellation can be expressed as [41], [42] 

^^e,MPSK = - / 07STBC (— ^^^j dO (147) 

TT Jq Vsm / 
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where = tt — vr/M and g = sin^ (vr/M). From (11471) . we can obtain the upper bound as 

Pe,MPSK< (l-j7]Krnci9;i) (148) 



which becomes tighter as 7 increases [49], and hence yields 

7^00 log 7 

Therefore, the asymptotic behavior of the MGF 0^5^.32 (s; 7) at large 7 reveals a high-SNR slope 
of the SEP curve. 

Suppose that 7 is sufficiently large. For nx < n^, it follows from (jV]) that 

log^^sTgc {g; 7) ^ - rank ^HJE!i$s ® ^t) ■ log7 + constant. (150) 

^ V ' 

?iT'™iii(nii,n.s) 

Similarly, using dH), we have for rix > n^, 

log0^sTBc (an) ~ - rank (^$r (g) HaH^) • log 7 + constant'. (151) 

^ V ' 

nH-min{nT,ns) 

Hence, 

(isTBC = min {m, n-n) ■ min {max (nx, wr) , ns} (152) 

from which dTUt follows immediately. For a general case of arbitrary two-dimensional signaling 
constellation with polygonal decision boundaries, the SEP can be written as a convex combination 
of terms akin to (11471) [82]. Hence, we can easily generalize the proof to the case of any two- 
dimensional signaling constellation. 

B. Proofs of Property \J\^ 

1) Proof of Property U} Let Ai,A2,...,A„ be the eigenvalues of $. Then, the correlation 
figure ( ($) defined in Definition |21 can be written as 

n 
k=l 

which is symmetric in Ai, A2, . . . , A„ and holds Schur's condition (1611) . Hence, we complete the 
proof. 

July 5, 2006 DRAFT 



42 REVISED FOR PUBLICATION IN THE IEEE TRANSACTIONS ON INFORMATION THEORY 

2) Proof of Property ^ Since Iliii C (^j) = C{®T=i^i)^ follows immediately from 
Property [TJ 

3) Proof of Property ^ Let A^*-*, A2'', • • • , A^*^ be the eigenvalues of $j (z = 1,2,..., m). 
Then, YlT=i C (^0 be written as 

i=l i=l k=l \ * / 

which is symmetric in |^A^*\ ^Ag^ • • • , ;]-An]}™^^ and holds Schur's condition (l6TT) . Since 
{^ll'^i^ ;II'^2^ • • • > ^-^"^^IHi the eigenvalues of we complete the proof. 

C. Proof of Theorem |21 

Using Theorem |7] in Appendix lUl we get 

= ' tr ($9 tr + tr ($9 tr (^^ + ' tr (^^ tr ($1) + {n^unf • 

(155) 

Combining (l24|), and (I155I) . together with the fact that E{||ir||p} = riTn-n, yields (ISTT) . 

D. Proof of Theorem 

In this case, the ergodic capacity (or Shannon-sense mean capacity) is given by the well-known 
expression [2]-[4] 

C (7) = E |log2 det (^In^ + ^HH^^ | bits/s/Hz (156) 

which is achieved by the complex Gaussian input X ~ J^nT,Nc{^nTxN^, -^IutiIn^)- 
From [53, (35)] and [53, Theorem 9], we get 



Eh_ ^ loge 2 ^ loge2 
A^omin E{||fr|||} r^R 



(157) 



DRAFT 



July 5, 2006 



SHIN AND WIN: MIMO DIVERSITY IN THE PRESENCE OF DOUBLE SCATTERING 



43 



and 



s 

E|tr ^(^fl"T)'J| 

2 (nrrn-onc)'^ 

(158) 



E=i =2 <; tr 



fet^^ 



Using Definition El and Theorem {7\ in Appendix HIl (11581 ) can be expressed in terms of the 
correlation figures of and $s as in ( B^ . 

Proof of Theorem |?| 

Due to the channel decoupling property of OSTBCs, the Shannon capacity of OSTBC MIMO 
channels can be written as 

CsTBC (7) = 7^ ■ E |log2 + ^^1^ j I bits/s/Hz (159) 

which is achieved by complex Gaussian inputs Xk ~ CA/'(0, ;^^^)- From [53, (35)], [53, Theo- 
rem 9] and the first two derivatives of (11591) at 7 = 0, it is easy to show (l44l) and 



27^ 

k{JH\ 



si^^ = ^,;;;„ ^ (160) 



If/ 

from which and Theorem |2l (l45t follows readily. 
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Fig. 1. Block diagram of a space-time block coded system in double-scattering (wt, ns, nR)-MIMO channels and induced 
SISO subchannels. 
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Fig. 2. SEP of 8-PSK (2.25 bits/s/Hz) versus 7 in spatially uncorrelated double-scattering (4, ns, 2) -MIMO channels, 
ns = 1, 2, 3, 5, 10, 20, 50, 100, cxd (i.i.d. Rayleigh). 
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SNR7 (dB) 

Fig. 3. SEP of 16-PSK Alamouti (4 bits/s/Hz) and Q4 (3 bits/s/Hz) OSTBCs versus 7 in spatially uncorrelated double-scattering 

(riT, ns, «.r)-M1MO channels. The Alamouti and Q4 codes achieve the diversity order of dsTBC = 2 in (2,3,1) and (4,2,1) 
links, respectively. The rfsxBc's for (2, 5, 3), (4, 3, 2) and (2, 10, 11), (4, 5, 5) pairs are 6 and 20, respectively. 
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SNR 7 (dB) 



Fig. 4. SEP of 8-PSK ^4 (2.25 bits/s/Hz) versus 7 in doubly correlated double-scattering (4, 10, 4)-MIMO channels. The 
transmit and receive correlations follow the constant correlation $t = *r = (p) for p = (spatially uncorrelated double- 
scattering), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9. For comparison, the SEP for i.i.d. Rayleigh-fading MIMO channels is 
also plotted. 
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Correlation coefficient p 



Fig. 5. SEP of 8-PSK Qi (2.25 bits/s/Hz) as a function of correlation coefficient p in doubly correlated double-scattering 
(4, ns, 4)-MIMO channels with constant correlation $t = *r = *4^' (p). '^s = 5, 10, 20, 50, 100, oo (doubly correlated 
Rayleigh) and 7 = 15 dB. 
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Fig. 6. SEP of 8-PSK Q4 (2.25 bits/s/Hz) versus ns in double-scattering (4, ns, 1)-MIM0 channels. The transmit and scatterer 
correlations follow the constant correlation *t = ^i"^ (p) and $s = *ns (p) for p = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 
and 0.9. 7 = 25 dB. 
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